Astronomy & Astrophysics manuscript no. Tycho 


©ESO 2013 


January 26, 2013 





Strong evidences of hadron acceleration in Tyclio's Supernova 

Remnant 

G. Morlino' * and D. Caprioli^'** 

INAF/Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125, Firenze, Italy 
Received August, 2011 ; accepted ~, - 

ABSTRACT 

Context. Very recent gamma-ray observations of G 120. 1 + 1.4 (Tycho's) supernova remnant (SNR) by Fermi-LAT and VERITAS 
provided new fundamental pieces of information for understanding particle acceleration and non-thermal emission in SNRs. 
Aims. We want to outline a coherent description of Tycho's properties in teims of SNR evolution, shock hydrodynamics and multi- 
wavelength emission by accounting for particle acceleration at the forward shock via first order Fermi mechanism. 
Methods. We adopt here a quick and reliable semi-analytical approach to non-linear diffusive shock acceleration which includes 
magnetic field amplification due to resonant streaming instability and the dynamical backreaction on the shock of both cosmic rays 
(CRs) and self-generated magnetic turbulence. 

Results. We find that Tycho's forward shock is accelerating protons up to at least 500 TeV, channelling into CRs about the 10 per cent 
of its kinetic energy. Moreover, the CR-induced streaming instability is consistent with all the observational evidences indicating a 
very efficient magnetic field amplification (up to ~ 300yuG). In such a strong magnetic field the velocity of the Alfven waves scattering 
CRs in the upstream is expected to be enhanced and to make accelerated particles feel an effective compression factor lower than 4, 
in turn leading to an energy spectrum steeper than the standard prediction oc E^^. This latter effect is crucial to explain the GeV- 
to-TeV gamma-ray spectrum as due to the decay of neutral pions produced in nuclear collisions between accelerated nuclei and the 
background gas. 

Conclusions. The self-consistency of such an hadronic scenario, along with the fact that the concurrent leptonic mechanism (inverse 
Compton scattering of relativistic electrons on several photon backgrounds) cannot reproduce both the shape and the normalization of 
the detected the gamma-ray emission, represents the first clear and direct radiative evidence that hadron acceleration occurs efficiently 
in young Galactic SNRs. 
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• 1. Introduction 

' The detection of high energy gamma-rays from supernova rem- 
I nants (SNRs ) has long been considered as the most promising 

■ tool (see e.g. iDrurv. Aharonian & Volk L Il994l) to probe the so- 
called supernova paradigm for the origin of Galactic cosmic rays 
(CRs), which states that SNRs are responsible for the accelera- 
tion ofnucldup to energies as high as a few times 10* GeV (see 

I e.g. lHmasll2005L for a review). 

In the last few years several SNRs have been detected in 
' gamma -rays, both in the GeV and in the TeV band (see Caprioli 
I ( 1201 lab for a recent summary) but, in spite of this larger and 

■ larger amount of data, a clear-cut evidence that such an emission 
is due to the decay of neutral pions produced in nuclear interac- 
tions between accelerated nuclei and the background plasma is 
still lacking. The very reason why the probe of such a hadronic 
scenario for the origin of the detected gamma rays is so contro- 
versial is that also processes involving relativistic electrons {lep- 
toni c scenario) may pro vide a similar gamma-ray signature (see 
e.g. Ellison e t alir2007l for a general discussion on this topic). 

The best example of this issue is indeed represented by the 
SNR RX J17 13.7-3946, which for many years has been re- 
garded as the best candidate as a hadron accelerator du e to its 
brigh t ness in the TeV band (see e.g. Morlino, Amato & Blasil 
l2009t IZirakashvili & Aharonianl 120101) . Nevertheless, very re- 
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cent high-re solution observatio ns in the GeV band performed by 
Fermi-LAT (lAbdo et al.L 1201 ih showed that the emission from 
the SNR shell is more likely due to Inverse-Compton scattering 
(ICS) of relativistic electrons rather than to pion decay. Its spec- 
trum in the GeV band is in fact very flat (oc E^'^-^), as predicted 
by a "leptonic" model, which is al so supported by the lack of 
detection of thermal X-ray emission lEllison et al. l (l20T0l) . 

Another leptonic mechanism which may be relevant 
for the gamma-ray emission from SNRs is the relativistic 
bremsstrahlung of accelerated electrons which, always coupled 
with ICS, has been proposed as a viable alternative to the 
hadronic scenario to explain for instance the spectrum of Cas A 
(Ab do et al.l 120101) . There are several other examples of SNRs 
whose emission could be, and has been, a ccounted for eith er in 
an hadronic or in a leptonic scenario (see ICaprioIii 1201 lal and 
references therein) but, at the moment, for no gamma-ray-bright 
source the leptonic scenario has been ruled out beyond a reason- 
able doubt. 

In several cases the extended TeV emission has been sug- 
gested to originate from the interaction between particles accel- 
erated in SNR and nearby molecular clouds (MCs), which natu- 
rally provide a large number of target nucleons. This scenario 
may be preferentially associated with remnants produced by 
core-collapse SNe, whose progenitors evolve very quickly in star 
formation regions still populated with parent MCs. In these en- 
vironments, however, it is not very clear whether to expect large 
CR acceleration efficiencies, because of the small ionized frac- 
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tion a nd, likely, because of ion-neutral damping (iDrurv et all 
Il996h . 

As pointed out for instance by ICapriolil (1201 lal) . the best 
sources where to look for hadron acceleration signatures are 
gamma-ray-bright SNRs not associated with MCs and showing 
rather steep spectral indexes both in the GeV and in the TeV 
bands. We identify Tycho's SNR, freshly detected by Fermi- 
LAT and VERITAS, as the best candidate for investigating the 
effects that particle acceleration produces at all observable wave- 
lengths, simultaneously using multiple set of data to constraint 
our model. To discuss things consistently, we nevertheless need 
to adopt a very refined model able to account, at the same time, 
for the SNR evolution, the acceleration of particles at the shock 
and, as it will be clearer soon, the magnetic field amplification 
produced by the streaming of the same relativistic particles. 

In this work we therefore apply the non-linear diffusive 
shock acceleration (NLDSA) theory to the Tych o's SNR follow- 
ing the semi - analytical approach put forward by lAmato & Blasil 
(l2005Ll2006h : ICaprioh. Amato & Blasil ( l2010blah and references 
therein. Thanks to the large amount of data available at vari- 
ous wave bands, this remnant can be considered one of the most 
promising object where to test the shock acceleration theory and 
hence the CR-SNR connection. Tycho is among the youngest 
Galactic SNRs, being only 439 years old: recent observations 
have in fact confirmed that it is the remnant of a Type la SN ex- 
ploded in 1572 by detecting its optical s pectrum near maxim um 
brightness from the scattered-light echo (iKrause et al 

Tycho shows a radio spectral index of 0.6-0.65 and a flux 
density as large as 40.5 Jy at 1.4 GHz (Kothes et al., 2006). Very 
interestingly, such a spectral index requires electrons in the GeV 
energy range to have a spectral slope ~ 2.2 - 2.3, i.e. steeper 
than the standard prediction of linear acceleration theory. The 
radio morphology is clearly shell-like, with enhanced emission 
along the northeastern edge (see Fig. [T]l. A hint of curvature in 
the radio spectrum, consistent w ith non-Unear theories' ex pecta- 
tions, has also been reported by Reynolds & EUisonl (Il992i) . who 
also inferred a magnetic field strength J 100//G. 

In gamma-rays, instead, Tycho has been observed several 
times with no credible detection reported until only very recently 
the VERITAS Cherenkov tel escope succeeded in measuring its 
TeV emission up to 10 TeV (lAcciari et al 1 1201 Ih . providing an 
integral flux over 1 TeV of only ~ 0.9% of the steady Crab 
Nebula emission above the same energy. Such a low emission 
makes Tycho one of the weakest TeV source ever detected. The 
TeV emission is compatible with a point like source since the 
telescope's point spread function (PSF) of ~ 15 arcmin is larger 
than the angular size of the remnant (8 arcmin). Even more re- 
cently, the Fermi Large Area Telescope (Fermi-LAT) has re- 
ported a detection of Tycho as well, assessing an integral flux 
above 400 MeV of ~ (3.5 ± + 0. 7,v,) x lO-^phcm-^s" ' 
and a spectral slope 2.3 ± 0.25,0/ ± 0. 1 r„r (lGiordano et aD. 1201 ll) . 

It has been raised the question whether the gamma-ray emis- 
sion were correlated with the interaction of the remnant with 
a nearby MC. The peak of the TeV emission might in fact be 
slightly displaced northeast to the remnant's center, where also 
a CO cloud has been detected, but the statistical significance 
of this displacement seems rather low (Acciari et al., 2011). A 
possible interaction between the northeastern part of Tycho and 
dense molecular gas has indeed been suggested by many au- 
fliors (iRevnoso et all 119991; iLee et all ^)04; Cai et al., 2009; 
IXu et alll20lT r but in a more recent studv lTian & Leahvi(i201 ll) 
analyzed the 21 cm continuum, the HI and '^CO-line data, con- 
cluding that there are no evidences of a direct interaction be- 
tween the shock and a MC. 



Interestingly enough, Tycho is also clearly visible in the 
infra-red (IR) band as well (llshihara et all l200ll) . In particular, 
mid-IR images show a shell-like structure with brightness peaks 
at the northeast and northwest boundaries, probably due to dust 
collisionally heated up by the hot shocked plasma in the SNR 
downstream. The far IR band, instead, is dominated by cold dust 
emission, which spatially correlates with the CO MC mentioned 
above, hence it is probably not related to the remnant. IR emis- 
sion due to heated dust is particularly important because it pro- 
vides the principal photons background for the ICS of acceler- 
ated electrons, as we will show in ^4.41 

High-resolution X-ray maps reveal strong non-thermal 



emi ssion con centrated in thin filaments (IHwang et al 



Bamba et al. 



2002 



201C; 



2005; W arren et all l2005t iKatsuda et al. 
lEriksen et all 1201 1 ). The thickness of these filaments is usu- 
ally interpreted as due to severe synchrotron losses of hig h en- 
ergy electrons radiating in a strong magnetic field (e.g. Bailed, 
l2006h . This interpretation allows us to estimate the magnetic 
field strength, which has to be as l arge a s ~ 250/iG inside the rim 
(IVolk et all 120051; iParizot et all |2006|) . There is however also 
the possibility that the observed X-ray morphology may result 
from a rapid dam p ing of the magnetic field in the downstream 
(iPohl et all bolOh . iMarcowifli & Cassel (l2010l) investigated this 
scenario considering different magnetic relaxation processes and 
different kind of magnetic turbulence and showed that, in order 
to have a relaxation length comparable to the X-ray rim thick- 
ness, the downstream magnetic field has invariably to be of order 
200 - 300;uG. Smaller fields would in fact require much longer 
distance to be damped. Hence, regardless of the presence of any 
damping mechanisms, the small thickness of X-ray filaments can 
be regarded as a proof that the magnetic field is amplified from 
the typical interstellar value up to hundreds of jiG. In ^2.41 we 
show that the presence of a damping mechanism for the mag- 
netic field can be investigated also by analyzing the radial profile 
of the radio emission. 

It is interesting to recall that detailed X-ray morphology 
studies have akeady provid ed evidences of effici ent hadron ac- 
celeration in Tycho. In fact IWarren et al.l (12005 ^ estimated the 
location of the contact discontinuity (CD) and the forward shock 
(FS) by comparing thermal and non-thermal X-ray emissions, 
concluding that they are too close to be described by purely 
gaseous hydrodynamical models. They therefore argued that the 
compression factor at the FS had to be larger than 4, as predicted 
in the case of efficient CR acceleration. 

Finally we note that hints on the presence of accelerated 
ions come also from the detection of Ho- lines. Analyzing some 
Balmer dominated filaments in the eastern limb of the remnant, 
iLee et al.l (1201 0) found a significant amount of the Ha emission 
to come from the region upstream of the shock, which has been 
interpreted as a CR-induced precursor 

All these things considered, several observations have al- 
ready provided strong hints that efficient CR acceleration is oc- 
curring at Tycho's FS. Our goal is therefore to provide a self- 
consistent picture able to take into account all of these observa- 
tional evidences and to couple them with brand new gamma-ray 
data. The NLDSA theory has already been adopted in order to 
predict the multi-wavelength spectrum of Tycho (e.g. lVoIk et all 
l2005l 120081) . nevertheless here we improve on previous attempts 
in some important aspects: 1) the magnetic field upstream of the 
shock is calculated from resonant streaming instability; 2) the 
dynamical reaction of the magnetic field onto the shock dynam- 
ics is accounted for; 3) the upstream scattering center velocity is 
properly calculated in the amplified magnetic field; 4) the non- 
linear Landau damping of the magnetic field in the downstream 
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is introduced; 5) the ICS of accelerated electrons is calculated 
considering as target photons non only the cosmic microwave 
background (CMB) radiation, but also the Galactic background 
and, more importantly, the IR photons produced by the local 
warm dust. 

The inclusion of the dynamical reaction of the field reduces 
the compressibility of the plasma and affects the prediction for 
the shock compression factor (Caprioli et al., 2009). A crucial 
ingredient is the velocity of the scattering centers, which is 
generally neglected with respect to the shock speed, but could 
be si gnificantly enhanced when th e ma gnetic field is ampli- 
fied ( Vladimirov, Ellison & Bykov, 2006' nCaprioli et al.L l2009t: 
Izirakasvh ili & Ptuskin. 2008)! When this occurs, the total com- 
pression factor felt by accelerated particles may be appreciably 
reduced and, in turn, the spectra of accelerated particles may be 
considerably softer. 

It is worth remembering that some observational features, 
especially the radio emission, are strongly affected by the past 
history of the remnant, hence any reliable calculation has to 
take into account also the SNR evolution. In this paper we 
use a stationary version of NLDSA theory, but we couple this 
theory to the hydrodynamica l evolution of the remnant pro- 
vided bv lTruelove&McKed ((1999) . We divide the SNR evo- 
lution in several time steps and we assume that for each time 
step the stationary theory can be applied, like has been done 
in [C aprioli. Amato & Blast (2010a). However, as showed by 
ICaprioli et al.l (I2OIOI) . stationary models and time-dependent ap- 
proaches return very similar CR spectra for non-relativist ic 
shocks. 

We compare the results of our kinetic model with the multi- 
wavelength integrated spectrum of Tycho from the radio to the 
TeV range, and also with the radial profile of X-ray and radio 
emissions. Our conclusion is that existing data of Tycho's SNR 
are consistent with a moderately efficient acceleration of CR nu- 
clei: at the present age we infer that a fraction around 12 per cent 
of the total kinetic energy has been converted in CRs. Such an 
efficiency also implies an amplified magnetic field of ~ 300;uG, 
perfectly consistent with the measured X-ray rim thickness. In 
addition, such a strong magnetic field enhances the velocity of 
the scattering centers, finally reducing the effective compression 
factor felt by accelerated particles, whose spectram turns out to 
be as steep as ~ E^^-^. The most important consequence of this 
fact is that this spectrum allows us to fit the observed gamma-ray 
emission, from the GeV to the TeV band, as due to neutral pion 
decay. Moreover, in this framework it is not possible to explain 
the TeV emission as due to ICS without violating many other 
observational constraints. 

The paper is organized as follows: in ^we summarize the 
details of our model for non-linear particle acceleration and our 
treatment of the SNR evolution. In ^ we outline the macro- 
scopic properties of Tycho's SNR, in order to fix the free param- 
eters of our model, while in ^we widely discuss the comparison 
between data and our findings for the multi-wavelength spec- 
trum, also by analyzing each different energy band separately. 
We conclude in ^ 

2. Description of the model 

2.1. Remnant evolution 

We model the evolution of Tycho by following the analytic pre- 
scriptions given by Truelove & Mc Kee ( 1999). More precisely, 
we consider a SN explosion energy Esn — 10^' erg and one 
solar mass in the ejecta, whose structure function is taken as 
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Fig. 1. Radio image of the Tycho's remnant at 1.5 GHz in lin- 
ear color scale. Image credit: NRAO/VLA Archive Survey, (c) 
2005-2007 AUI/NRAO. 
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Fig. 2. Time evolution of shock radius Rsh, shock velocity Vsh, 
magnetic field immediately behind the shock B2 and CR accel- 
eration efficiency = Pcr/puV^/,- 



(y/vej) ^ (see §3.2 and §9 in lTruelove & Mc Ked.[T99l . Such 
a set of parameters has been showed to be suitable for describ- 
ing the evolution of the FS position and velo city for a type la 
SNR: the parametrization given in table 7 of Trti elove & Mc Keg 
(Il999h in fact differs from the exact numerical solution of about 
3 per cent typically, and of 7 per cent at most. Such a solution, 
which does not include explicitly the possible role of the CR 
pressure in the SNR evolution, is still expect to hold for moder- 
ately small acceleration efficiencies (below about 10 per cent). 
We checked a posteriori that the efficiency needed to fit obser- 
vations does not require amore complex treatment of the shock 
evolution during the ejecta-dominated stage. 

The circumstellar medium is taken as homogeneous with 
proton number density no = 0.3 c m'-^ and temperature T o = 
10"* K. Following the conclusion of lTian&Leahvl dloTlh . we 
assume that the remnant expands into the uniform interstellar 
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medium (ISM) without interacting with any MC. With these pa- 
rameters, the reference value for the beginning of the Sedov- 
Taylor stage is Tst = 463yr, therefore Tycho is experiencing 
the transition between the ejecta-dominated and the adiabatic 
stage. FS radius and veloc ity can be hence simply written as 
dTruelove & Mc Keel [19991) : 
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The time evolution of some relevant physical quantities is 
depicted in Fig.|2] The last point coiTesponds to an age of 439 yr 
and hence to our estimate for the present shock velocity, Vsi, - 
4990 km s"', and shock radius, Z?,/, = 3.94 pc. 

The radial structure of density and temperature profiles is 
then calculated by assuming that the downstream (from the CD 
to the FS), being subsonic, is roughly in pressure equilibrium. In 
our case, this recipe leads to a discrepancy of less than 10 per 
cent with respect to exact profile calculated by including also 
the contribution of t he CR pressure in the SNR evolution (see 
e.g. lChevaheil[T98l . 

More precisely, since p''(f) cc pgasit) cc V^i^(t), the adiabatic 
decompression of a fluid element with density po shocked at time 
to can be calculated as 



po 



V,h{t) 



Vshito) 



(3) 



where y is the effective adiabatic index of the CR+plasma fluid. 
In general we have 4/3 < y < 5/3, but in our case y = 5/3 is a 
very reasonable choice, justified by the fact that the CR acceler- 
ation efficiency we find is never much larger than ~ 10 per cent. 
The quantity L(fo, 0^1 introduced above represents the adia- 
batic energy loss CRs and magnetic fields undergo in the period 
t - to, and it is consistently taken into account in modelling the 
emission from fluid elements shocked at earlier times. 

In a similar way, this simple recipe allows us to track also 
the downstream temperature profile and therefore to work out 
the bremsstrahlung emission due to thermal electrons, which in 
general have a temperature T,, different from the protons' one, 
Tp. We assume = Tp nie/mp immediately after the shock and 
calculate the heating of electrons downstream as due to Coulomb 
collisions against hotter protons. In Fig. |3]we plot the value of 
electron and proton temperatures between the CD and the FS at 
the present epoch, normalized to the proton temperature imme- 
diately behind the shock, Tp ^h = 5 x 10** K. The figure shows 
that Coulomb collisions rapidly enhance the electron tempera- 
ture up to few per cent of Tp ,,^, and also that the maximum value 
is reached close to the CD, where we have ^ 2 x 10^ K. We 
would like to stress that Coulomb collisions provide the mini- 
mum heating rate for electrons; others plasma processes could 
in principle equilibrate protons and electrons even more rapidly. 
On the other hand, as we will see in ^4.21 the thermal electron 
bremsstrahlung predicted by assuming Coulomb heating only is 
compatible with the non de tection of X-ray thermal emission 
dCassam-Chenai et al.Lr2007h . 

2.2. Particle acceleration 

On top of this SNR evolution, the spectrum of accelerated par- 
ticles is calculated acco rding to the semi-analyti c kinetic for- 
malism put forward in ICaprioli. Amato & Blasil (1201 Obi) and 
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Fig. 3. Density (solid line) and temperature (dashed line) of the 
shocked ISM in the region between the CD (vertical solid line) 
and the FS, computed at the present age of the remnant, t - 439 
yr, and normalized to the values immediately behind the shock. 
The dot-dashed line shows the electron temperature computed 
including the heating due to Coulomb collisions. 



references therein (especially lAmato & Blasil (l2005l l2006l) X 
which solves self-consistently the equations for conservation 
of mass, momentum and energy along with the diffusion- 
convection equation describing the transport of non-thermal par- 
ticles for quasi-parallel, non-relativistic shocks. In particular, 
we impose the CR distribution function to vanish at a dis- 
tance ~ XescRsh upstream of the shock, mimicking the presence 
of a free-escape boundary beyond which highest-energy parti- 
cles cannot diffuse back at the s hock and get lost in the ISM 
dCaprioh. Amato & Blasil 12010b'). This constraint actually de- 
termines the instantaneous maximum momentum than can be 
achieved by CRs by imposing the diffusion length of particles 
with momentum Pmaxit) to be comparable with the distance be- 
tween the shock and the free escape boundary, namely 



y, 



sh 



" Xescl^sh 



(4) 



where Dip) is the difiiision coefficient, here taken as Bohm-like, 
i.e. D(p) - ^-^riip), with v{p) and riip) the particle velocity 
and the Larmor radius, respectively. 

We also checked a posteriori that by posing Xesc -O.l such a 
size-limited pmax is also consistent with requiring both the accel- 
eration time up to Pmax to be smaller than the age of the system 
(iBlasi. Amato & Caprioli. 2007) and the diffusion length down- 
stream to be smaller than the distance between FS and CD. 

Particles are injected into the acceleration mechanism from 
the suprathermal tail of the Maxwell ian distribution (thermal 
leakag e) following the prescription by iBlasi. Gabici & Vannonil 
(l2005h . We assume that all the particles with a momentum 
larger than ^inj times the downstream thermal momentum have a 
Larmor radius large enough to cross the shock from downstream 
to upstream, and in turn to be injected in the acceleration pro- 
cess. We fix ^inj - 3.7, corresponding to let a fraction ~ lO"'* of 
the particles crossing the shock to be accelerated as CRs. 

It is worth stressing that the solution of the NLDSA prob- 
lem obtained with the semi-analytical approximate solution 
worked out in Caprioli, Amato & Blasi (2010b) is in excellent 
agreement with the exact solution of the stationary diffusion- 
convection equation for CR transport coupled with conservation 
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Fig. 4. Total spectrum of protons and electrons (divided by K^,,) 
at the present age of the remnant, as a function of the particle 
momentum. For comparison, the dot-dashed line shows a power- 
law oc p "^^ which emphasizes a small flattening of the proton 
spectrum at the highest momenta. 



equations, but it is also very consistent with state-of-the-art nu- 
merical approaches to the time-dependent problem and also with 
Monte Carlo approaches able to retain all the information about 
the anisotropy of the whole (therm al particles + CRs) distribu- 
tion function (Caprioh et al.ll2010h . 



2.3. Magnetic field amplification 

The super-Alfvenic streaming of relativistic particles is ex- 
pected to excite a plasma instability (streaming instability) 
which strongly enhances Alfv en waves resonant wi t h the L armor 
radii of relativistic particles (ISkillingl Il975t iBelll Il978h . Such 
a magnetic turbulence may grow well beyond the quasi-linear 
limit, eventually producing a perturbation 6B » Bq, where Bq 
is the component of the ordered background magnetic field par- 
allel to the CR gradient, and hence to the shock normal. This 
process, usually referred to as magnetic field amplification, has 
rather solid observational counterparts in the X-ray observations 
of many young SNRs: the downstream narrow X-ray rims due 
to severe synchrotron losses by relativistic electrons (see e.g. 
IVolk et all 12001 iParizot et al.[ 12006), the ra pid variability of 
the X -ray-bright spots in RX J1713. 7-3946 dUchivama et all 
l2007h and a lso the limited exten sion of the upstream emission 
in SN 1006 (iMorlino et al.Ll2010l) . 

Tycho makes no exception, showing evidences of narrow 
non-thermal X-ray rims pointing to magn e tic fields as l arge as 
a few hundreds lu G (Hwa ng et al.L 12 002: Bambaetal., 2005; 
IWarren et al.ll2005; Katsuda et al.Ll2010.:.Eriksen et al., 2011) . 

We model magnetic field amplification d ue to resonant 
streaming instability as in ICaprioli et al ] (120091) . i.e., by assum- 
ing that saturation is achieved when P„ ^ Pcr/2MA, where P„. 
and Per are the pressure in Alfven waves and in CRs, normalized 
to the ram pressure poV^f^, and Ma is the Alfvenic Mach number 

Ma - VshlvA - Vsh yj47Tpo/B. In the Hmit My » 1,M^ » 1 
the solution for the wave transport equation reads (see eq. 43 in 
ICaprioli etaLll2009l) : 



1 + Ujx) 
AMa{x)U(x) 



Pcr(x), 



(5) 



where U(x) is the fluid velocity (normalized to Vsi,) in the up- 
stream precursor 

When magnetic field amplification occurs a strong non- 
linear regime, as it is expected in our case, it is not clear 
whether the magnetic field entering the Alfven velocity has to 
be the background one, Bo, as predicted by an extrapolation 
of the quasi-linear theory, or the amplified one, 6B » Bo- 
Strictly speaking, the quasi-linear prediction is related to the fact 
that resonant transverse Alfven waves (6B ± Bq) of arbitrary 
strength are exact solutions of the problem: it is however very 
likely that some mechanism may tend to make the field rather 
turbulent when 6B » Bo, therefore we choose to calculate the 
Alfven velocity in the amplified magnetic fi eld. The implicat ions 
of this choice are discussed more widely in ICaprioTil (|20 1 1 ab . 

Here and in the following we label with the subscript 1 
(2) quantities measured immediately ahead (behind) the shock, 
while the subscript is reserved for undisturbed quantities at up- 
stream infinity. The maximum value of the upstream magnetic 
field is thus reached immediately ahead of the shock and it is 
given by 



Bi 



1 -I- U 



Vl',1 



(6) 



where we posed ^cr = Pcr.i - I - Ui. The relevant Alfvenic 
Mach number thus reads 



M, 



A.i 



2 {l-^crY 



l-U f ^cr 2- ^„ 

and therefore we have upstream of the subshock: 



B 



Ma 



(l-^c.F^ 



(7) 



(8) 



where we also used the mass conservation pi Ui - po- 

We then assume the downstream magnetic field to have a 

strength B2 = ^(1 -h 2/?^ J/3 Bj ^ VTTBi, where R,„h ^ 4 



is the subshock compression ratio, since only the components 
perpendicular to the shock normal are compressed. For the 
shock parameters considered above (po - 0.3mpCm"-', V,/, ^ 
5000 km s ') a CR acceleration efliciency of about 5-10 per cent 
is therefore enough to provide a downstream magnetic field of a 
few hundreds yuG (see Eq. [8]l, in agreement with the values in- 
ferred by X-ray observations of young remnants. 

As outlined in the introduction, such large magnetic fields 
have additional consequences on both the shock hydrodynam- 
ics and on the spectrum of accelerated particles. More precisely, 
when the upstream pressure in magnetic fields becomes com- 
parable to, or even larger than, the thermal plasma pressure, 
the global dynamics of the system is regulated by the interplay 
betwe en the CR pre s sure and the pressure in magnetic turbu- 
lence (ICaprioli et al 1 12OO8II2OO9I) . To be more quantitative, at 
the shock Pvv.i - 1/2M2 ^ ^ 2.3 X IQ-^ for^„ = 0.06 and hence 
Ma.i - 15, which is much largerthanfgj ^ l/jM^ ^ 5.6x10"^ 
since the sonic Mach number is Mo - 443. The additional pres- 
sure and energy density in magnetic field thus help to balance 
the compression of the upstream fluid induced by CR pressure 
and is consistently included both in the conservation equations 
and in the cal culation of the jump condi tions at the subshock, as 
described bv Caprioh et all (12008112009 1). Moreover, such a mag- 
netic feedback helps in keeping the compression ratios rather 
close to 4, which in turn legitimates the adoption of Eqs. (Tjand 

m 
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The second, and probably most important, effect of the effi- 
cient magnetic field amplification is that the phase velocity of the 
Alfven waves CRs scatter against may become a non-negligible 
fraction of the fluid velocity. In particular, Alfven waves pro- 
duced by the CR gradient via streaming instability travel in op- 
posite direction with respect to the fluid, so that upstream CRs 
actually feel a net velocity u = u - va - u{l - 1/Ma), while 
downstream it is likely for helicity mixing to lead to v^o x 0. 
Since for DSA the spectral slope is determined only by the com- 
pression ratio f ^ (ui - va.i)/u2 = mi/m2(1 - ^/Ma) felt by 
diffusing particles, the non-negligible velocity of the scattering 
centers (often also called Alfvenic drift) naturally leads to par- 
ticle spectra steeper than what predicted without including this 
effect. 

The potential role of Alfvenic drift in the DSA theory has 
already been pointed out by Bell (1978), but it is important to 
stress that steep particle spectra can be produced at SNR shocks 
only if the magnetic field is amplified via some inst ability to 
levels corresponding to Ma - 20 or lower (see also ICapriolil 
l20TTal) . In our case, taking the saturation of the streaming insta- 
bility in the amplified field makes this phenomenon even more 
non-linear, with the final result that the larger the number of 
particles injected, the stronger the amplified magnetic field and 
eventually the steep er the CR spectr a achieved in a consistent 
NLDSA calculation (CaprioliL lIoTTbl) . 

Finally, another fundamental aspect of magnetic field ampli- 
fication is that diffusion in the enhanced magnetic turbulence al- 
lows particles to return to the shock more quickly, in turn achiev- 
ing a much larger momentum both in the age-limited and in the 
space-limited scenario (see ^2. 2b . In our calculations we there- 
fore take the Bohm-like diffusion coefficient in the amplified 
magnetic field Bi rather than in Bo, finding p^n-^^ as high as ~ 500 
TeV for the accelerated protons. 

Nevertheless, the excitation of modes resonant with CR gy- 
roradii is not the only phenomenon which may be re s ponsi- 
ble for mag netic field amplification. As sho wn by IBelll (|2004|) 
and then by iRiquelme & Spitkovskvl (l2010l) . also non-resonant 
modes with short wavelengths may effectively grow because 
of CR escaping the system. The relative importance of reso- 
nant and non-reson ant modes has been inv estigated within a ki- 
netic approach by lAmato & Blasil (l2009l) . who found that for 
young, strong shocks also non-resonant modes may play a non- 
negligible role. A comprehensive treatment of the interplay be- 
tween CRs and magnetic fields from a kinetic point of view is 
however beyond the goal of this paper, therefore we simply as- 
sume a uniform magnetic field in the upstream given by Eq. dS). 
Such a recipe is indeed very reasonable in the shock proxim- 
ity and becomes more and more heuristic in the far upstream, 
where amplification should be prevalently due to the excitation 
of Bell's modes. However, if far upstream the magnetic field 
were much lower than we assumed, we would have two effects: 
on one hand, the CR spectrum would be slightly flatter (oc £ '^^) 
at the highest energies, where particles would feel a total com- 
pression ratio of about R,ot - 4.2 and, on the other hand, the 
maximum energy would be consistently reduced. These two ef- 
fects would indeed have a radiative signature in the TeV band 
but the quite large measurement errors and the lack of detection 
of a high-energy cut-off do not allow us to make a realistic in- 
vestigation of a possible flattening at the highest energies, yet. 

In addition, also long-wavelength modes may play a role in 
the problem enhancing the diffusion of the highest-energy parti- 
cles (see, e.g., Bvkov et al., 2011, and references therein) up to 
PeV fractions. Our simple choice, very common in the literature, 
of adopting a constant magnetic field upstream is indeed an over- 



simplification of the problem but, within our incomplete under- 
standing of the mechanisms responsible for the scattering of the 
highest-energy particles, allows us to heuristically account for 
a crucial observational evidence: if the GeV-to-TeV emission is 
hadronic, as we will demonstrate in the following, protons have 
to be accelerated up to at least a few hundreds TeV (see also 
lEriksen et al.Ll201 ll) . 

2.4. Magnetic field damping 

As it will be clearer in the following discussion, in order to repro- 
duce both the intensity and the radial profile of the synchrotron 
emission, it is necessary to account for some magnetic field 
damping in the downstream as well. Several damping mech- 
anisms have been proposed to be effective in SNRs (see e.g. 
Ptuskin & Zirakasvhili, 2003, and references therein), but here 
we focus only on non-linear Landau damping, which is expected 
to be rat her efficient in hot (low Mach number) plasmas, as o ut- 
hned bv lVolk & McKenzid lll981h : lMcKenzie & VoUd <1982h . 

We therefore assume that in the downstream the magnetic 
field is damped with a rat e, averaged on the different scales, 
given by (see eqs. 10-12 in lPtuskin & Zirakasvhiliil2003l) : 



r„, - 0.05 ^„i„VA -0.05 



VA ^ 0.05 va(Bi)c 

rdPmax) 3 XescVshRsh 



(9) 



where we took the maximum turbulence scale (i.e. the minimum 
wave number k^m) to be given by the Larmor radius of particles 
with momentum close to p^ax (in the last equality we simply 
made use of Eq. lU. 

Since B\Ui - B2U2, the typical length-scale for the non- 
linear Landau damping is thus given by 



U2Vsh 3Xesc ^l, 



0.05 VAiB2)c 



Rsh * 3.5 PC 



(10) 



and therefore we obtain the following recipe for the downstream 
magnetic field: 



B(r) - B2 exp 



R ,h - r 



(11) 



The actual radial profile of Bif) is finally showed in Fig. |5] 
(solid line) and compared with the undamped magnetic field 
(dashed line) and with the magnetic field subjected only to adia- 
batic expansion (dotted line). 

2.5. Electron spectrum 

Once the proton spectrum has been calculated self-consistently 
with the shock dynamics, it is also possible to work out the spec- 
trum of accelerated electrons, which is parallel to the protons' 
one since DSA is charge independent, but truncated at a lower 
maximum energy determined by synchrotron losses in the am- 
plified magnetic field. 

The instanta neous electr on spectrum at the shock, fe,oip), is 
computed as in lMorlino . "Am ato & Blasi (200 91), using th e ex- 
pression worked out by IZirakashvili & AharonianI (|2007|) for a 
oc proton spectrum: 



fe.oip) = Kep fpfi(p) 



1 -H 0.523 (p//?,,„ax)^ 



-p-/pI, 



(12) 



where pe.max is determined by equating the acceleration time 
with the energy loss timescale due to synchrotron emission. 
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Electrons produced at time to with e nergy £0 will thu s have 
an energy E(t) at a later time t given bv (lRevnoldsLll998l) : 



0.92 0.94 

r/Rsi, 



1.00 



Fig. 5. Radial profile of the magnetic field in the downstream at 
the present age of the remnant. The dashed curve shows B2 as 
result from the amplification and compression at the shock. The 
dotted curve include adiabatic losses, while the solid curve is 
computed including both the adiabatic expansion and the Landau 
damping. 



which, as a consequence of the efficient magnetic field amplifica- 
tion, turns out to be much small er than Tycho's age for large mo- 
menta (see Eqs. (3) and (4) in iMorlino. Amato & Blasj 120091) . 
Such a spectrum is expected to be a reasonable approximation 
of the true one (seee.g. lBlasill2010h . smce our proton spectrum 
is almost a power-law only slightly steeper than p"^, as showed 
in Fig. m 

A very important property of Eq. (ITZi is also to retain the 
correct shape of the cut-ofF, which is not a simple exponential, as 
assumed by s everal authors, but rather goes as exp [-(/?/ Pe max)^] 
jBlasiL 12010) . This trend is in fact crucial when studying the X- 
ray emission in the keV band. 

Finally, the constant Kep accounts for the different normal- 
ization between electron and proton spectra, very likely related 
to the different mechanisms responsible for lepton and hadron 
injection. However, since electron injection in SNRs is still far 
from being u nderstood from first principles (though see e.g. 
lAmano & Ho shino. 2010; Riquelme & Spitkovsky, 2010) we do 
not fix Kgp a priori but we leave it as a free parameter which has 
to be tuned in order to fit the observations, and in particular the 
synchrotron emission. 

Another important piece of information needed to describe 
the non-thermal emission due to relativistic electrons is their 
downstream evolution. After being accelerated at the shock, 
electrons are in fact advected downstream, losing energy be- 
cause of the remnant adiabatic expansion and of radiative losses 
(mostly synchrotron, but also ICS). The standard equation de- 
scribing the energy evolution from to to t in this case reads; 



dE 4 
at 3 



Stt 



m 

Eo 



L{to,t) 



1 +AEo J' L{to,T)Bl/T)dT^ 



(14) 



EdL 
Z d7' 



(13) 



where A = 1.57 x lO"-' in cgs units. 

The electron spectrum at time t can be therefore computed 
by using number conservation, namely: fe(E)dE - fe o{Eo)dEo. 
The spectra of accelerated protons and electrons integrated over 
the downstream volume, at the present age of the remnant, are 
finally showed in Fig. |4] and in particular the electron spectrum 
can be described as follows. 

For p i proii - 200 GeV/c, where is the roll-over mo- 
mentum, i.e. the momentum for which the radiative loss time 
becomes comparable with the SNR age, synchrotron losses are 
negligible and so electron and proton spectra are parallel and 
oc p"^-^. Between proii and pcui - 20 TeV/c, i.e. the momentum 
for which the loss rate exceeds the acceleration one, the electron 
spectrum is instead dominated by radiative losses, and therefore 
goes as ~ p^^-^. 

It is interesting to notice that the steepening due to syn- 
chrotron losses spans more than one order of magnitude in mo- 
mentum below pcui, therefore it is expected to have peculiar ob- 
servational counterparts both in the synchrotron and in the ICS 
emissions. We will discuss this topic in details in ^4.3l and 34.4! 
respectively. 

2.6. Radiative processes 

In order to fit the emission observed from Tycho, we consider 
the following processes: 1) synchrotron emission of relativistic 
electrons; 2) thermal and non-thermal electron bremsstrahlung; 
3) ICS of electrons on different microwave, IR and optical pho- 
tons; 4) photons due to the decay of neutral pions produced in 
hadronic collisions. 

The synchrotron emission produced by accelerated electrons 
in the local amplified magnetic field (Eq. [TTI) is carried out by 
using the exact synchrotron kernel (Rybicki & Lightman, 1985|). 

The bremsstrahlung emission consists of a component due 
to thermal electrons, whose temperature is calculated as in ^2.51 
plus a non-thermal component produced by the accelerated ones. 
The total spectrum is the sum of both electron-nucleus and 
electron-electron bremsstrahlung: latter process is however im- 
portant at the highest energies, being negligible at lower ones. 
For the e- e and the e-n bre msstrahlung w e adop t the differen - 
tial cross sections provided bv lHaud (1 19981) and bv lHaugl([T997h . 
respectively. 

Electrons contribute to the gamma-ray emission through 
ICS on local photons. We consider here the CMB radiation, 
the IR dust emission and the Galactic IR + optical back- 
ground as a function o f the distance from the Galactic center 
jPorter & Strong! |2005|) . The flux and spectrum of ICS photons 
is calculated by using the exact kernel for ICS, with also the full 
Klein-Nishina regime accounted for. 

Finally, flux and spectrum of the gamma rays produced by 
7r° decay are calculate d following the analytical approximations 
bv lKelneretan(l2006h . 



where L - L(to, t) accounts for adiabatic losses, as in Eq. (O, and 



B^-hB£„ withBe 



^gq, ^eq - 12yuG is the equivalent magnetic field 
which takes into account the IC losses due to CMB radiation. 
Galactic background Ught and IR emission due to local dust (see 



3. Modelling Tycho 

In this section and in the next we apply the whole apparatus out- 
lined above to explain Tycho's properties, along with its multi- 
wavelength spectrum and morphology as well. 
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Our model actually has very few free parameters other than 
the well-constrained ones related with a type la SN explosion: 
the number density of the upstream medium, no - Pa I nip, the 
injection efficiency ^inj and the electron to proton normalization. 
Kg p. The values of no and ^inj are simultaneously chosen in or- 
der to fit the gamma-ray emission detected by Fermi-LAT and 
VERITAS, in such a way that no also allows us to reproduce 
the inferred position and velocity of the SNR forward shock. 
Eventually, we may regard also the choice of;^'^,^. = 0.1 and the 
assumption of a constant magnetic field equal to B\ upstream as 
a heuristic way to account for the maximum hadron energy and, 
in turn, for the highest-energy photons detected by VERITAS 
(see the discussion at the end of ^2. 3b . 

The best fitting is obtained by adopting no - 0.3 cm"^ 
and finj = 3.7. The chosen value for the upstream density 
is compatible with the upp er limits existing in the literature: 
ICassam-Chenai' et 2A\ (l2007h derived no < 0.6 cm from the 
absence of thermal emission from the shocked ejecta in the 
Chandra X-ray data, while a sim ilar result, no < 0.4 cm"'', 
was obtained by ' Volk et al] (1200 8') using the upper limits on the 
gamma-ray emission measured by HEGRA and Whipple. 

The value of no can be checked indirectly also by compar- 
ing the current estimates of the remnant's distance and of the 
shock velocity with the same quantities computed according to 
the evolution model explained in § |2] Being Tycho the remnant 
produced by a type la SN, we can safely fix the explosion en- 
ergy in EsN - 10^' erg and the ejecta mass in about one solar 
mass. Therefore, once also the density of the interstellar medium 
is fixed, all other quantities related to the SNR evolution are un- 
equivocally determined. 

Our estimate for the distance can be worked out by compar- 
ing the final radius of the FS with the observed size of the rem- 
nant: using no = 0.3 cm"^ we get d - 3.3 kpc. In the literature 
different techniques have been adopted to provide an estimate of 
such a dis tance, always return ing a value ranging between 2 and 
5 kpc (see lHavato et al.Ll20101, for a review). Kinema tic methods 
give a distance of 2.5 - 3 kpc (iTian & Leah vll20 111) . It is worth 
noticing that kinematic estimates are affected by the fact that 
Tycho is located in the Perseus arm of the outer Galaxy, which 
is influenced by the spiral shock which causes a velocity rever- 
sal. Other methods tends to prefe r higher distan ces. In particular 
using the detection of echo light. iKrause et alj (2008) determine 
a distance of 3.8^Qg kpc. Combining the observed ejecta veloc- 
ities with the ejecta proper-motion measurements by Chandra, 
Havato et alj (2010) derived a distance of 4 ± 1 kpc. Finally 



Volk et al.l (12008.) suggested a lower limit of 3.3 kpc by mod- 



eling the existing gamma-ray measurements. In summary our 
result of 3.3 kpc seem to be consistent with all the existing esti- 
mates. 

The shock speed that we found at the current age of the rem- 
nant is Vsh - 4990 km s Measuring the proper motion of 
X-ray rim observed with the RO SAT saXeWiit, Hughes (2010) de- 
rived Vsh - 4600 ± 400((i/2.3 kpc) km s"', in perfect agreement 
with our findings. Similar results have been ob tained measur- 
ing th e proper motion in the radio frequencies (iRevnoso et al.L 
11997 *). Another independent technique sometimes used to infer 
the shock speed is the measurement of the optical Ha emission. 
In fact the line width of the so called broad Ha emission is re- 
lated to the temperature of the shocked plasma, w hich is in turn 
related to the shock speed ( Chevalier et al.L [T980I) . By adopting 
this technique Smith et al.l (1 19911) derived shock velocities in the 
range 1500-2800 km s"', about a factor 2-3 lower than our es- 
timate. Nevertheless, we have to notice that the shock speed in- 
ferred from the Wa line width is systematically lower than the 



ones provided by other methods in basically all the SNRs consid- 
ered. This discrepancy is probably due to the fact that the down- 
stream temperature is calculated using a regular hydrodynamical 
shock, without accou nting for the possibl e presence of acceler- 
ated particles (see e.g.l Helder at al.Ll2009l) and/or the dynamical 
role of the neutral component. 

Finally, we fix the electron to proton ratio from the fit of the 
synchrotron emission. Very interestingly, a unique value Kep - 
1 .6x10"^ allows us to fit both the X-ray dat a from Suzaku and the 
radio emission from several experiments (iRevnolds & EUisonL 
1992). Our value of K^p is only a factor a few smaller than 
the ratio measured at Earth in the diffuse spectrum of Galactic 
CRs, and such a discrepancy might be accounted for by remem- 
bering that the latter value is the result of all the SNR evolu- 
tion (in principle K^p is a function of time) and it is also av- 
eraged over the contributions of different kinds of sources (e.g. 
type I/II SNRs and maybe pulsars as well). At the present time 
there are some hints a bout how electron injection into PSA 
may occur dAm ano & Hoshinoi 120101: iRiquelme & Spitkovskvl 
2010': 'M orUnoi 1200 9), but a complete theory coupled to NLDSA 
and therefore to a self-consistent calculation of the non-thermal 
emission from a SNR during its lifetime is still lacking. This 
is actually the very reason why two of our free parameters are 
somewhat related to the need of modeling the injection of parti- 
cles into the acceleration mechanism. 

With both model and environmental parameters fixed, we 
find that the Tycho's FS is converting into CRs about the 6 per 
cent of its bulk pressure (Fig. |2]i, which in turn implies a down- 
stream magnetic field B2 - 300yL(G. Our shock dynamics is only 
mildly modified by the presence of accelerated particles: apart 
from magnetic field amplification, all of the signature typical of 
NLDSA, like the formation of an upstream precursor where the 
incoming fluid is slowed down (and in turn the concavity in the 
spectrum due to the fact that CRs with different momenta probe 
different regions of the precursor) are not very marked. For in- 
stance, the total and the subshock compression ratios turn out to 
be R,oi - I/U2 - 4.2, and Rsub = U1/U2 - 3.9, respectively. 
This result also justifies a posteriori the treatment of the SNR 
evolution as governed by the gas pressure only 02. lb rather than 
by the whole gas -1- CRs fluid. 

Let us now compare the theoretical prediction for the ratio 
between the positions of C D and FS with the val ue inferred from 
the X-ray observations by Warr en et al.l (l2005l ) which returned 
RcolRsh = 0.93 ± 2%. In our model we find Rcd/R^i, = 0.87. 
This value has however to be taken with a grain of salt, since 
it may be enhanced by Rayleigh-Taylor instabilities at the CD. 
For instance, purely gaseous ID hydrodynamical simulations 
predict a significantly small er value for this ratio, and namely 
RcolRsh = 0.77 ( Wang & Chevaheil [2001). In the same work 
the authors also performed a 2D simulation and showed that the 
Rayleigh-Taylor instability at the CD allows fingers of ejecta to 
protrude well beyond the average CD radius, also inferring the 
maximum extent of these fingers in about the 87 per cent of the 
FS radius. Following this line of thought, we can estimate in 
about a 5 per cent the boost in the CD/FS position ratio due to 
instabilities operating at the CD. When we multiply our predic- 
tion by 1.05, we get a CD/ FS ratio of a bout 0.91, consistently 
with the value estimated W Warren et al.l(l2005l). 

Some authors (e.g. lEUison. Decourchelle & BalletL l2005l) 
have proposed that young SNRs can accelerate CRs efficiently 
also at the reverse shock (RS). This idea is supported by the 
fact that the RS moves with respect to the unshocked ejecta 
with a velocity comparable to that of the FS. The RS speed 
can in fact be estimated from the usual analytical model by 
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Fig. 7. Surface brightness of the radio emission at 1.5 GHz as a 
function of the radius (data as in Fig.[TJ. The thin solid line rep- 
resents the projected radial profile computed from our model, 
while the thick solid line corresponds to the same profile convo- 
luted with a Gaussian with a PSF of 15 arcsec. 



iTruelove & Mc Ke3 (Il999l) . and for Tycho it turns out to be 
~ 0.86 Vsh- On the other hand, an efficient DSA requires also 
a sizable magnetic field, say of order of few //G, to efficiently 
scatter particles. The simple dilution by flux-freezing of the typ- 
ical magnetic field on the surface of the progenitor would imply 
~nG or lower magnetic field, therefore a very efficient magnetic 
field amplification at the RS would be requir ed to make DSA 
work. On the other hand, IWarren et al.l (l2005l) showed that the 
inferred distance between the CD and RS is compatible with hy- 
drodynamical model excluding efficient CR acceleration at the 
RS. 

Anyway, in case of non negligible acceleration at the RS, we 
would expect an enhancement of the synchrotron emission only 
from the position of the RS towa rds the inside of th e remnant. 
According to the measurement of IWarren et al.l (l2005l) the RS is 
located at OJRsh- From Fig.|2]we do not notice any clear excess 
in this region with respect to our prediction. We will discuss in 
^4.11 how also our analysis is able to account for the observed 
radio profile without requiring acceleration of electrons at the 
RS. 



4. Multi-wavelength spectrum 

In Fig. |6] we show our best fitting of the photon spectrum pro- 
duced by the superposition of all the radiative processes outlined 
above, comparing it with the existing data. The overall agree- 
ment is quite good, therefore we want to analyze now in greater 
detail the emission in each single band. 

4.1. Radio emission 

As depicted in Fig.|6l the total radio emission is nicely accounted 
for by our model, which returns a photon spectral index of .606 

in the considered energy range (10 to 1500 MHz). 

Ne vertheless, as already noticed by iRevnolds & Elliso"nl 
(|1992|) . Tycho's radio emission shows evidence for a curvature 
in the spectram, which turns out to be slightly steeper (flatter) 
than E-^-^ below (above) 100 MHz. In fact, Kothes et al. ( 2006) 
find a spectral index a = 0.65 for the radio spectrum fitted to 
all flux densities and a = 0.61 in the range 408 to 1420 MHz. 



IRevnolds & ElUsonl (Il992h ascribed this effect to the curvature 
of the electron spectrum due to non-linear effects induced by the 
CR precursor and stressed the fact that, in the energy range in- 
volved in the radio emission, electron and proton energy spectra 
are not exactly parallel because protons are non-relativistic and 
DSA is momentum dependent. 

In our model we did not find any steepening in the lowest 
energy region of the electron spectram, and the reason is the 
following. The emission in the 1 - 100 MHz band, where the 
curvature is observed, is produced by electrons with energy close 
to the injection energy. The typical synchrotron frequency is in 
fact: 



AeBy^ 



= 24 



B 



100/iG 102 



-^MHz, 



(15) 



from which, taking the typical injection Lorentz factor ye,mi ~ 
10 - 100 and a downstream magnetic field strength of ~ 100 - 
300yL(G, we obtain Vsyn ~ 1-100 MHz. The shape of the electron 
spectrum in this energy range is therefore quite uncertain, and 
our simple parametrization of particle injection via thermal leak- 
age may indeed be too simplistic. A comparison between ther- 
mal leakage model and more refined Monte Carlo approaches 
to particle injection, able to retain the smooth transition in the 
supra-thermal energy range, seems however to confirm that a 
proper account for particle injection should quite naturally pre- 
dict a spectral steepening in the 10-100 MeV energy range (see 
fig. 1 of CaprioH et al., 2010). 

Therefore, the detection of a curv ature in the low-ener gy ra- 
dio emission from young SNRs (Revnolds & EllisonL[T99"2 ) may 
not be necessary related to the presence of a CR precursor in the 
upstream, but rather be a ciTicial tool for probing electron injec- 
tion. It is also worth stressing that this possibility is given by 
the fact that the magnetic field is strongly amplified: for a stan- 
dard interstellar field, in fact, the same radio emission would 
come from higher-energy electrons, and more precisely in the 
GeV range. 

Besides the volume-integrated emission, another precious in- 
formation can be infeiTed from the remnant morphology: the ra- 
dial profile of both radio and X-ray emissions provides in fact 
a strong evidence that the magnetic field inside the SNR is con- 
siderably amplified. In order to compare the observed radial pro- 
files, we need to compute the emission projected along the line 
of sight. As already noticed, Tycho is clearly shell-like both in 
radio and X-ray bands: even if the northeastern region shows 
an enhancement of the emission, especially in the radio band, 
the spherical symmetry can be indeed considered a reasonable 
approximation. Under this assumption the local emissivity is a 
function of the radius only and the projected emission is thus 
simply due to the integration along the line of sight /: 



jpiy,p) = 2 



(16) 



the integration limits take into account that the emission is ex- 
pected to come only from the shocked ISM located between the 
CD and the FS. 

In order to be compared with the actual data, the projected 
emission has to be convolved with the instrumental point spread 
function, which we assume to be Gaussian. We use the radio 
map in Fig.[T] which has been obtained with the VLA on January 
2007 at the frequency of 1.5 GHz, (data from the NRAO/VLA 
Archive) and in Fig. Q we compare our results with these data. 
The step line shows the measured surface brightness obtained 
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Fig. 6. Spatially integrated spectral energy distribution of Tycho. The curves show synchrotron emission, thermal electron 
bremsstrahl ung and pion decay as ca lculated within our model (see text for details). The experimental data are, respectivley: 
radio from jR evnolds & EHisonl (1 19921) : X-rays from Suzaku ( courtesy of Toru T amagawa) , GeV gamma-rays from Fermi-LAT 
dGiordano et al.. 2011.) and TeV gamma-rays from VERITAS (lAcciari et all l20i 11 Both Fermi-LAT and VERITAS data include 
only statistical error at 1 cr. 



by integrating the emission along the azimuthal angle, between 
and In. The thin solid line shows the projected emission com- 
puted using Eq. (fTSl l. while the thick line is the emission after the 
convolution with the instrumental point spread function, which 
is equal to 15 arcsec. 

Fig. |7] shows indeed a good overall agreement between the 
data and our prediction, even if some little discrepancies can 
be noticed. The most evident one is that in the innermost re- 
gion {r/Rsi, i 0.3), the theoretical prediction overestimates the 
observed brightness by about 20 per cent. A plausible explana- 
tion of this difference may reside in a slight deviation from the 
spherical symmetry, which is somehow expected just because 
the northeastern region is brighter than the rest of the remnant. 

Another subtle but interesting difference is that the emis- 
sion peaks slightly more inwards than in our model; as a con- 
sequence, also the emission detected in the region 0.6 i r/Rsi, ~ 
0.8 is found to be a bit larger than the theoretical prediction. 
This difference might have different explanations. The most ob- 
vious, and already mentioned, is the possible deviation from the 
spherical symmetry. Another possibility is given by placing the 
CD in a different position: if one assumed the CD to be located 
closer to the center (i.e. if one took the CD/FS ratio to be a few 
per cent smaller), the theoretical prediction would nicely fit the 
data. However, we can not forget that this explanation would be 
at odds with the findings of Warren e t al. (2005) . who estimated 
the position of the CD to be more towards the forward shock, 
namely around 0.93/? 



A final comment on the radio profile concerns the effects of 
the non-linear Landau damping in the determination of the mag- 
netic field relevant for the synchrotron emission. If we neglected 
the damping, the magnetic field strength in the downstream (dot- 
ted line in Fig.|5]l would lead to a total radio flux larger by a fac- 
tor 50 per cent or more with respect to the data, even if the radial 
radio profile would retain a rather similar shape. 

4.2. X-ray emission 

As it is clear from Fig.|6] the synchrotron emission spans from 
the radio to the X-ray band, where it sums up with the emission 
due to thermal bremsstrahlung. 

The best-fitting to the X-ray continuum observed by Suzaku 
data is illustrated in greater detail in Fig. |8] where the dashed line 
indicates the synchrotron emission alone and the solid line cor- 
responds to the sum of synchrotron plus thermal bremsstrahlung. 

The electron temperature in the downstream, calculated tak- 
ing into account only the heating due to Coulomb collisions with 
protons (Fig. [3]l, results in a bremsstrahlung emission peaked 
around 1.2 keV which, at its maximum, contributes for about 
the 6 per cent of the total X- ray continuum emission onl y, in 
agreement with the findings of ICassam-Chenai' et al.l (l2007h . 

In the same energy range there is however a non-negligible 
contribution from several emission lines, which becomes more 
and more important moving inwards from the FS, where the X- 
ray emission is mainly non-thermal (IWarren et al.l l2005l) . A de- 
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Fig. 8. X-ray emission due to synchrotron (dashed line) and to 
synchrotron plus thermal bremsstrahlung (solid line). Data from 
the Suzaku telescope (courtesy of Toru Tamagawa). 
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Fig. 9. Projected X-ray emission at 1 keV. The Chandra data 
points are from Cassam-Chenai' et al. (2007) (see their Fig. 15). 
The solid line shows the projected radial profile of synchrotron 
emission convolved with the Chandra point spread function (as- 
sumed to be 0.5 arcsec). 



tailed model of the line forest is, however, beyond the main goal 
of this paper. 

The projected X-ray emission profile, computed at 1 keV, is 
shown in Fig. [9l where it is compared with the Chandra data in 
the region that lCassam-Chenai et all (l2007l) call region W. The 
solid curve represents the resulting radial profile, already con- 
voluted with the Chandra PSF of about 0.5 arcsec, and shows a 
remarkable agreement with the data. As widely stated above, the 
sharp decrease of the emission behind the FS is due to the rapid 
synchrotron losses of the electrons in a magnetic field as large 
as ~ 300;uG. In Fig.|9]we also plot the radial radio profile com- 
puted without magnetic damping (dashed line); since the typical 
damping length-scale is ~ 3 pc, it is clear that the non-linear 
Landau damping can not contribute to the determination of the 
filament thickness. 

It is worth stressing that the actual amplitude of the mag- 
netic field we adopt is not determined to fit the X-ray rim profile, 
but it is rather a secondary output, due to our modelling of the 
streaming instability, of our tuning the injection efficiency and 



12 14 
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Fig. 10. Synchrotron emission calculated by assuming constant 
downstream magnetic field equal to 100 (dotted line), 200 
(dashed line) and 300 /iG (solid line). The normalization of the 
electron spectrum is taken to be Kej, = 1.6 x 10"^ for all the 
curves. 



the ISM density in order to fit the observed gamma-ray emis- 
sion (see the discussion in |3]l. We in fact checked a posteriori 
whether the corresponding profile of the synchrotron emission 
(which, in shape, is also independent on Kep), were able to ac- 
count for the thickness of the X-ray rims and for the radio profile 
as well. 

4.3. Radio to X-ray fitting as a liint of magnetic field 
amplification 

Another very interesting property of the synchrotron emission is 
that a simultaneous fit of both radio and X-ray data may provide 
a downstream magnetic field estimate independent of the one 
deduced by the rims' thickness. 

In fact, assuming Bohm diffusion, the position of the cut-off 
frequency observed in the X-ray band turns out to be indepen- 
dent of the magnetic field strength, actually depending on the 
shock velocity only. 

On the other hand, if the magnetic field is large enough to 
make synchrotron losses dominate on ICS and adiabatic ones, 
the total X-ray flux in the cut-off region depends only on the 
electron density, in turn fixing the value of Kep independently 
of the magnetic field strength. Moreover, radio data suggest the 
slope of the electron spectrum to be equal to 2.2 at low energies, 
namely below £■,„// ^ 200 GeV. Above this energy the spectral 
slope has in fact to be 3.2 up to the cut-off determined by setting 
the acceleration time equal to the loss time, as discussed in ^2.51 

In Fig.[TO]we plot the synchrotron emission from the down- 
stream, assuming a given magnetic field at the shock and ne- 
glecting all the effects induced by damping and adiabatic expan- 
sion. The three curves correspond to different values of B2 = 
100, 200 and 300;uG, while the normalization factor K^p is cho- 
sen by fitting the X-ray cut-off and it is therefore the same for all 
curves. As it is clear from the figure, in order to fit the radio data 
the magnetic field at the shock has to be J 200/iG, even in the 
most optimistic hypothesis of absence of any damping mecha- 
nism acting in the downstream. 

As a matter of fact, synchrotron emission alone can provide 
an evidence of ongoing magnetic field amplification, indepen- 
dently of any other evidence related to X-ray rims' thickness or 
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emission variability. Such an analysis is in principle viable for 
any SNR detected in the non-thermal X-rays for which it is also 
possible to infer the spectral slope of the electron spectrum from 
the radio data, only requiring radio and X-ray emissions to come 
from the same volume and therefore from the same population 
of electrons. 

4.4. Gamma-ray emission 

The most intriguing aspect of Tycho's broadband spectrum is its 
gamma-ray emission, which has been detected before in the TeV 
band by VERITAS (Acciari et al., 20111) and then in the GeV 



band by Fermi-LAT, too (.Giordano et al.L 1201 ih . 

Gamma-ray emission from SNRs has been considered for 
long time a possible evidence of hadron acceleration in this class 
of objects (Drurv, Aharonian & Volk , 1994), even if there are 
two distinct physical mechanisms which may be responsible for 
such an emission: in the so-called hadwnic scenario the gamma- 
rays are produced by the decay of neutral pions produced in nu- 
clear collisions between CRs and the background gas, while in 
the so-called leptonic scenario the emission is due to ICS or rel- 
ativistic bremsstrahlung of relativistic electrons. 

We show here, with unprecedented clarity for a SNR, that the 
gamma-ray emission detected from Tycho cannot be of leptonic 
origin, but has to be due to accelerated hadrons, instead. This 
fact, along with the VERITAS detection of ~ 10 TeV photons 
and the lack of evidence of a cut-ofF in the spectrum, implies 
that hadrons have to be accelerated up to energies as high as a 
few hundreds TeV. 

In particular, the proton spectrum we obtain shows a cut- 
ofF around /?max - 470 TeV/c (see Fig. HJi. In this respect, 
Tycho could be considered as a half-PeVatron at least, be- 
ing there no evidence of a cut-ofF in VERITAS data. The 
age-old problem of detecting SNRs emitting photons with en- 
ergies larger than a few hundreds TeV (i.e., responsible for 
the acceleration of particles up to the knee observed in the 
spectrum of diffuse Galactic CRs) may therefore be only a 
matter of time. Of course it may also be that actually not 
all the SNRs are PeVatrons, or even that it is quite unlikely 
to observe a SNR in the brie f stage of its life in which it 
may be acting like a Pe Vatron (^ Caprioli. Blasi & Amatol 120091: 
ICaprioli. Amato & Bla si. 2010a). 

This very important result relies on the fact that the spectral 
slope in the sub-GeV to multi-TeV range is consistent, within 
measurement errors, with a single power- law oc ii"^ ' - E^^-^. 
Such a power-law is steeper than the test-particle prediction 
of DSA and, more importantly, is significantly different from 
what expected in the framework of standard NLDSA theories, 
which predict conca ve spectra fl atter than E^^ a t the high- 
est energies (see e.g. Drurvl 19831: Blandford & Ei chleH, 119871: 



Malkov & Drurvi i200h for some re- 



iJones & ElhsonL 11991 
views). 

A key role in our calculations is in fact played by hav- 
ing included in the NLDSA also magnetic field amplification, 
and in particular by having assumed the relevant Alfven veloc- 
ity as the one calculated in the amplified magnetic field. The 
fact that a small Alfvenic Mach number may imply quite steep 
spectra had already been put forward by iBelll (11978 ). but only 
by coupling this fact with a very efficient streaming instability 
allows NLDSA theories to predict a proton spectrum steeper 
than E^^ over the whole energy range (Zirakasvhili & Ptuskinl 
I2OO8; C aprioli et al .. 2009), with also the possibility of chang- 
ing the spectral slope by tuning the injection efficiency dCaprioUl 
1201 Ihh . 



In this respect, fixing a value for ^inj regulates both the slope 
and the normalization of the proton spectrum and, in turn, of the 
gamma-ray emission. In calculating the photon spectrum pro- 
duced by :/!■" decay (Fig.fTTI) we also accounted for the presence 
of He nuclei in the target medium and for the likely presence of 
accelerated particles heavier than protons, which are expected to 
boost the total emission by a factor ~ 1 .4 and ~ 1.5, respectively 
(Mori, 2009; Caprioli, Blasi & Amata l20HT) . 

There are two main reasons why the observed gamma-rays 
cannot be of leptonic origin: 1) the intensities of both ICS and 
relativistic bremsstrahlung emission are too low and 2) the ex- 
pected shape of the leptonic emission is radically different from 
the one observed. 

The former point is illustrated in Fig. [TT] where we com- 
pare gamma emissions due to bremsstrahlung, ICS and pion de- 
cay. The contribution due to relativistic bremsstrahlung (dot-dot- 
dashed line), however, is negligible at all the wave-lengths of 
interest. The ICS has been computed for three different target 
photon families: CMB radiation. Galactic background radiation 
and IR due to local dust emission. The ICS on the CMB radiation 
(dashed line) provides a contribution a factor 20-30 smaller than 
the one by pion decay (dot-dashed line), while the ICS on the IR 
-I- optical Galactic background (dotted line) at Tycho's position 
is even smaller. F or the Galactic background we used the esti- 
mate provided by Porter & Strong (2005) for a distance from the 
Galactic Center of 12 kpc and inside the Galactic plane, which is 
a good approximation for Tycho's position (~10 kpc away from 
the Galactic center). 

For the ICS we have also included, for the first time as far as 
we know, a third contribution due to local dust emission (solid 
thin line). In fact IR radiation from Tycho's direction has been 
detected by different satel htes in the waveleng th band ranging 
from 10 jim to 160 //m ("I shihara et al.L 12001 1) . This emission 
seems to be due to two distinct populations of dust, one which 
we refer to as warm, with a temperature T ~ 1 00 K and one 
we ref er to as cold, with T ~ 20K. According to llshihara et al.l 
(1200 ih . the warmer component should originate from ISM dust 
collisionally heated by the hotter plasma downstream of the 
shock. In fact, the emission in the mid-IR band shows a shell- 
like structure recalling both synchrotron and Ha emissions. On 
the other hand, the far-IR emission (140-160 fj.m) is dominated 
by a rather cold dust emission correlating with the position of the 
CO cloud, which has been showed not to be in physical contact 
with the SNR itself. Hence we consider only the photon back- 
ground due to the warm dust. In order to compute the local en- 
ergy density, we assume that IR photons are emitted uniformly 
from the shocked I SM between the CD and the FS. Following 
llshihara et al.l (12001 ) we fix the temperature of warm dust to be 
100 K, hence, by fitting the observed IR data, we get a local 
energy density of Ewarm =3.1 eV/cnr'. It is important to notice 
that the ICS on this mid-IR radiation dominates over all other 
photon backgrounds. Furthermore we notice that if one assumed 
also the far-IR emission (due to cold dust) to originate from a re- 
gion close to the remnant, the ICS on these latter photons would 
be negligible with respect to the ICS associated with warm dust 
photons. Summing up all the possible photon backgrounds, the 
resulting total ICS emission contributes not more than ~ 20 per 
cent of the total gamma emission around 0. 1 TeV, but it is negli- 
gible for lower and higher energies, as clearly shown in Fig. [TT] 

The possibility of explaining the gamma emission using only 
the ICS on the CMB photons by assuming a larger value of Kep 
is inconsistent with the very efficient magnetic field amplifica- 
tion inferred. A K^p smaller by a factor 20-30 would in fact re- 
quire a downstream magnetic field a factor 5-6 smaller in order 



G. Morlino and D. Caprioli: Strong evidences of hadron acceleration in Tycho's Supernova Remnant 



13 



VERITAS 




Total 



10^ 



10^ 



10 



10 



10 



11 



10 



12 



10 



13 



E [eV] 



Fig. 11. Gamma-ray emission observed by Fermi-LAT and by VERITAS compared with spectral energy distribution produced by 
pion decay (dot-dashed line), relativistic bremsstrahlung (dot-dot-dashed) and ICS computed for three different photon fields: CMB 
(dashed). Galactic background (dotted) and IR photons produced by local warm dust (solid). The thick solid line is the sum of all the 
contributions. Both Fermi-L AT and VERITAS data points include only statistical errors at Icr. For VERITAS data the systematic 
error is found to be ~ 30% (lAcciari et all 1201 L) . while for Fermi-LAT the systematic uncertainties are comparable or even larger 
than the statistical err or especially for the lowest energy bins due to difficulties in evaluating the galactic background (see Fig. 3 in 
IGiordano et al.Ll2bl ll and the related discussion). 



to keep the synchrotron emission fixed, in turn raising issues not 
only with the radial profile of the synchrotron emission (see also 
^4.2b . but also with the relevant effects induced by a large Alfven 
velocity discussed in 32.31 

The second solid argument against a leptonic origin for the 
detected gamma-ray emission is given by the actual slope of the 
photon spectrum. Once excluded the underdominant relativis- 
tic bremsstrahlung and ICS on the CMB radiation and Galactic 
background, we are left with ICS on the IR background due to 
local dust as the only viable candidate. However, as predicted 
by standard ICS theory and as showed in Fig. [TT] the expected 
photon spectrum below the cut-off is typically flatter than par- 
ent electrons' one, and more precisely is oc v"" for an electron 
, clearly at odds with Fermi-LAT data in the 



-2.2 



spectrum oc E 
GeV range. 

Another point worth noticing is that the ICS on the CMB 
radiation is sensitive to the steepening of the total electron spec- 
trum above ~100 GeV (Fig. |4| due to the synchrotron losses 
particles undergo while being advected downstream, while for 
the ICS on the IR+optical background the onset of the Klein- 
Nishina regime (above » 7 TeV for photons of 1 eV) does 



not allow us to probe significantly the steep region of the elec- 
tron spectrum. 

In other words, ICS on the CMB radiation is too low and 
cannot be boosted by invoking a larger electron density, while 
ICS on IR and/or optical background, which might as well be 
locally enhanced with respect to the mean Galactic value, cannot 
provide a spectral slope in agreement with both Fermi-LAT and 
VERITAS data. 

We are therefore forced to conclude that the present multi- 
wavelength analysis of Tycho's emission represents the best ev- 
idence of the fact that SNRs do accelerate protons, at least up to 
energies of about 500 TeV. The proton acceleration efficiency is 
found to be ~ 0.06poV^/,, corresponding to converting in CRs 
a fraction of about 12 per cent of the kinetic energy density 
^PoV^^i^. As estimated for instance in §3 of the review by Hillas| 
( 2005), such a value is consistent with the hypothesis that SNRs 
are the sources of Galactic CRs, provided that the residence time 
in the Milky Way scales with ~ E^^^^. 

It is important to remember that the actual CRs produced 
by a single SNR is given by the convolution over time of dif- 
ferent contributions with non trivial spectra, and namely the 
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flux of particle escaping the remnant from upstream during 
the Sedov-Taylor stages and the bulk of particles released in 
the ISM at the SNR 's d eath ( CaprioH, Blasi & Amato, 2009; 
ICaprioli. Amato & BlaisiL l2010al) . In this respect, the instanta- 
neous spectrum of accelerated particles in Tycho, which is in- 
ferred to be as steep as oc E^^-^, provides a hint of the fact that 
SNRs can indeed produce rather steep CR spectra as required 
to account for th e oc E^^-'' diffuse spectrum of Galactic CRs 
(ICaprioliLl20TTbl) . 



5. Discussion and conclusions 

Thanks to the large amount of data available at different wave- 
lengths, Tycho's SNR is one of the best object where study- 
ing the connection between CRs and SNRs. Using the observed 
non-thermal spectrum, and in particular the recent detection of 
GeV emission by Fermi-LAT together with the TeV spectrum 
detected by VERITAS, we can infer that at Tycho's FS protons 
are accelerated up to energies as large as ~ 500 TeV, and that the 
total energy converted into CRs can be estimated to be about 12 
per cent of the FS bulk kinetic energy. 

To reach this conclusion we investigated particle accelera- 
tion at the forward shock using a state-of-the-art semi-analytical 
NLDSA model including the dynamical reaction of the accel- 
erated particles, the generation of magnetic fields as due to 
streaming instability excited by CRs, the dynamical reaction of 
these self-generated fields on the plasma and also the modifi- 
cation of the speed of the scattering cente rs (Alfvenic drift) in- 
duced by the magnetic field amplification (ICaprioU et all l2009t 
[Caprioli, Amato & Blasi, 2010a). The last efffect is of crucial im- 
portance because it produces a softening of the particle spectrum 
with respect to the standard prediction oc E'^, which allow us to 
fit both the radio and the gamma photon spectra. The stationary 
version of NLDSA we used is coupled to the dynamical evolu- 
tion of the remnant according to the analytical prescription by 
iTruelove & Mc Ked (fl999h . This allow us to account for infor- 
mation sensitive to the time evolution, like for instance the radio 
to X-ray connection, as illustrated in ^4.31 

We would like to stress that our fit to the observed non- 
thermal spectrum and to the SNR morphology depends on a 
few parameters only. As for modeling particle acceleration, we 
are forced to introduce some parameters describing the non- yet- 
understood processes related to the microphysics of the magnetic 
fields involved, like the injection of hadrons and electrons (two 
parameters, ^,„y and K^p) and the scattering of highest-energy 
particles {xesc, plus a heuristic choice of a constant magnetic 
field upstream, B\). 

Being Tycho the remnant of a type la SN explosion, the only 
free parameter related to the SNR environment is actually the 
density of the circumstellar medium where the remnant expands 
(see We adopt a number density of no = 0.3 cm"^, which is 
consistent with all of the observational constraints coming from 
thermal and non-thermal emission and which also provides po- 
sition and velocity of the FS in decent agreement with the esti- 
mates provided by several authors adopting different techniques. 

Our findings show that the decay of neutral pions produced 
in hadronic collision between accelerated ions and gas nuclei is 
the dominant process in the gamma-ray band. In particular, we 
predict a slope for accelerated protons oc E^^^ which well ac- 
counts for Fermi-LAT and VERITAS detections within the ex- 
perimental errors. The ICS of relativistic electrons cannot ex- 
plain the observed gamma-ray emission, as can be seen from 
Fig.[TT]and as explained in 34.41 



The very reasons supporting this conclusion are the fol- 
lowing. First, the strong magnetic field produced by the CR- 
induced streaming instability forces the number density of rel- 
ativistic electrons to be too small (the electron to proton ratio is 
Kep = 1.6 X lO"-') to explain the gamma-ray emission as due to 
ICS on the ambient photons. A larger Kgp would in fact lead to 
overestimate the synchrotron emission, both in the radio and in 
the X-ray bands. 

Second, even if we arbitrarily reduced the magnetic field 
strength, enhancing at the same time the electron number density 
in order to fit the TeV gamma-rays with ICS emission, we could 
not account for the GeV gamma-rays because both the spectral 
slope and the flux would be incompatible with the recent Fermi- 
LAT observations. 

Also the other competing leptonic process which may show 
the correct spectral slope throughout the whole gamma-ray spec- 
trum, namely the non-thermal bremsstrahlung, has to be ruled 
out, in that it provides a flux two order of magnitudes lower than 
the Fermi-LAT detection, and cannot be arbitrarily enhanced 
without overpredicting both the TeV and the X-ray emission. 

In this work we also showed, for the first time as far as we 
know, that the main contribution to the ICS is due to the IR pho- 
tons produced by circumstellar dust heated up through collisions 
with the shocked plasma, which dominates on the CMB and 
Galactic photon background. This contribution, which is gener- 
ally neglected in multi-wavelength studies, may be relevant also 
for other SNRs, especially those produced by type la SNe which 
expand in the cold ISM. 

As we outlined above, the impossibility of fitting the gamma- 
ray spectrum with leptonic processes is tightly connected with 
the fact that many observational evidences point towards a very 
efficient magnetic field amplification, as we comment in the fol- 
lowing. 

In our model CR-induced streaming instability amplifies 
the magnetic field upstream of the shock from the unperturbed 
Galactic value of ~ 5^G up to ~ 90//G, which becomes ~ 300/iG 
immediately downstream because of compression at the shock. 
Observationally speaking, such a large magnetic field has two 
main consequences: 1) it produces narrow X-ray filaments be- 
cause electrons undergo severe synchrotron losses; 2) it deter- 
mines the roll-over frequency of the synchrotron spectrum to be 
~ 4 eV/fi, (see 34.31 1. The former effect has been extensively in- 
vestigated in the literature and it is indeed recovered in our study 
of the X-ray emission profile (see Fig. |9]l. 

Also the latter effect represents an independent evidence that 
the magnetic field has to be amplified up to a few hundreds of 
jiG. As we showed in § 14.31 when the radio data are good enough 
to infer the electron spectral index, the value of Vr„u, i.e. the fre- 
quency emitted by electrons whose loss time equals the age of 
the remnant, is uniquely determined by a simultaneous fit of both 
radio and X-ray data as due to synchrotron emission. In fact, 
when synchrotron losses are dominant, v™// is a function of the 
average magnetic field downstream {B) only, and from Fig. [10] 
we can see that for Tycho we have Vrou ~ 4 eV/h, which in turn 
implies {8} ^ IQQfiG. 

The narTowness of X-ray filaments might be produced also 
by rapid damping of magnetic field downstrea m of the shock 
rather than by synchrotron losses (see e.g. Poh l et al.l l2010h . In 
the case of Tycho we can exclude this possibility simply by look- 
ing at the synchrotron emission: if the field were damped below 
200 ;uG on a length-scale comparable to the X-ray rim width, 
we would in fact fail to fit the combination of radio plus X-ray 
emission, and the radio emission profile as well. This conclusion 
of course does not imply any damping not to be at work in the 
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downstream, but only suggests that no damping mechanisms can 
determine the structure of the observed X-ray filaments. 

Similar structures in Tycho's interior have been recently de - 
tected also in the shape of radial stripes by Eri ksen et al.N201 ih . 
These stripes show a surprisingly ordered pattern and have been 
interpreted by the authors as evidences of magnetic fields of few 
hundreds fiG, compatible with acceleration of particles up to the 
knee. Nevertheless, we would like to stress that a large mag- 
netic field is necessary but not sufficient to claim protons to actu- 
ally be accelerated up to consistently large energies: we showed 
NLDSA at the FS to be able to account for protons of about 
500 TeV, i.e. with an energy only a factor a few lower than the 
knee's in the CR spectrum. Larger energies might in principle 
be achieved at earlier or later times in Tycho, or even in peculiar 
region of the FS with a favorable magnetic field topology. In any 
case, we cannot exclude either the possibility that type la SNRs, 
like Tycho, might not be responsible for the knee obsei-ved in 
the diffuse spectrum of Galactic CRs, or the possibility that we 
are somehow underestimating E,„a_^, since VERITAS data do not 
provide evidence for a high-energy cut-off. 

Some authors have proposed that CR acceler- 
ation may occur also at the RS of yo ung SNRs 



(see 



Ellison. Decourchelle & Balletl 



2005 



Zrrakashvili & Aharoniani 2010t IZirakasvhil i & Ptuskin 



201 Ih . even if it is unlikely to expect a very large magnetic field 
in the ejecta RS propagate into. In the present work we therefore 
investigated the radial profile of the radio emission, concluding 
that there are no evidences of acceleration at the RS. In ^4.11 
we showed that it is possible to account for the observed radio 
profile including electrons accelerated at the FS only. 

In this work we tried to use basically all of the available 
observations, interpreting them in the light of a semi-analytical 
approach to NLDSA to build up a consistent model of particle 
acceleration at Tycho's FS. However, we are aware that we ne- 
glected the information coming from the detection of Balmer 
lines, which also point towards efficient particle acceleration 
(iLee et all l20Toh . The reason is that, at the moment, a reliable 
NLDSA theory in presence of neutrals still lacks. When shocks 
propagate into a partially ionized circumstellar medium, as it 
might be the case for Tycho, the neutral fraction cannot be ne- 
glected because it is coupled to the ions via charge-exchange and 
ionization processes. As a consequence, shock dynamics and, in 
turn, particle acceleration are likely sensitive to this aspect of the 
problem. The study of these phenomena might indeed provide 
new insights on how particle acceleration works, but, neverthe- 
less, it is unlikely for these effects to invalidate the most impor- 
tant findings of the present work, namely that Tycho is provid- 
ing us with the first clear-cut example of hadron acceleration in 
SNRs, up to ~ 500Ge V at least. This result, made possible by the 
present generation of gamma-ray instruments as Fermi-LAT and 
VERITAS, is of primary importance in corroborating the SNR 
paradigm for the origin of Galactic CRs. 
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